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摘 要 :和 裂纹 识别 是 结构 健康 监测 的 重要 内 容 。 本 研究 基于 反 演 分 析 原 理 ,将 数值 流 形 方法 
(NMM) 5 Elman 神经 网 络 相 结合 开展 裂纹 识别 。NMM 用 于 获取 对 应 裂纹 构 型 下 测 点 的 位 移 数据 
以 供 Elman 神经 网 络 的 学 习 , 在 此 基础 上 利用 训练 好 的 Elman 网 络 进行 了 直线 裂纹 反 演 。 通 过 2 
个 典型 算 例 证 实 了 NMM-Elman 协同 方法 的 可 行 性 和 精度 ,与 此 同时 分 析 了 测 点 布置 方式 及 输入 数 
据 噪 声 等 因素 对 裂纹 反 演 精度 的 影响 。 表 明 本 研究 的 方法 能 够 准确 反 演 出 单一 及 复杂 裂纹 的 裂 央 
间 玲 标 。 本 研究 的 工作 为 复杂 裂纹 的 高 效 准确 识别 提供 了 一 种 新 的 思路 和 方法 。 
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Crack inverse analysis with the NMM-Elman neural 
network collaborative method 
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'Abstract ; Crack identification is an important issue in structural health monitoring. Based on the principle 
“of inverse analysis, this paper combines the numerical manifold method ( NMM) with the Elman neural 
network to carry out crack identification. To serve the learning of Elman neural network, the NMM is used 
to obtain the displacement data of measuring points under various crack configurations. On this basis, the 
trained Elman network is used for straight crack inversion. The feasibility and accuracy of NMM-Elman 
collaborative method are verified by two typical examples. At the same time ,the effects of measuring point 
layout and input data noise on crack inversion accuracy are analyzed. The research shows that the method 
proposed in this paper can accurately reflect the crack tip coordinates of single and complex cracks. This 
work provides a new pathway for efficient and accurate detection of complex cracks. 
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裂纹 是 结构 体 中 常见 的 缺陷 ,其 存在 会 对 结构 。” 红外 热 成 像 技 术 、 压 电阻 抗 技 术 、CT 技术 等 ) ,但 其 
的 可 靠 性 产生 不 同 程度 的 影响 。 在 过 去 几 十 年 中 ， ”或 多 或 少 存在 检测 效率 低 、 成 本 高 .适用 面 罕 、 精 度 


已 经 发 展 了 许多 裂纹 无 损 检测 技术 ( 如 
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馈 声 波 技术 、 不 佳 等 问题 ;此 外 ,上 述 方法 一 般 需 预 佑 裂纹 的 大 体 
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位 置 ,导致 测 区 有 限 而 常 出 现 漏 检 的 情况 , 留 下 工程 
隐患 。 因 此 ,发 展 和 研究 新 型 无 损 检测 方法 对 结构 
体 的 健康 监测 和 寿命 维护 具有 重要 意义 。 

近年 来 ,数值 计算 工具 和 机 器 学 习 算法 的 飞速 
发 展 为 缺陷 检测 提供 了 另 一 种 路 径 。 此 类 方法 一 般 
先 通过 数值 仿真 技术 ( 如 有 限 元 法 、 扩 展 有 限 元 法 、 
等 几何 方法 .比例 边界 元 法 等 ) 进行 正 分 析 以 获得 结 
构 体 的 受 荷 响应 ,在 此 基础 上 利用 各 类 智能 优化 算 
法 (如 人 工 神经 网 络 .遗传 算法 、 粒 子 群 优化 Jaya 算 
法 等 ) 对 缺陷 进行 反 演 。 文 献 [1] 提 出 了 一 种 基于 
红外 热 波 技术 、 有 限 元 与 支持 向 量 机 相 结 合 的 复合 
材料 缺陷 检测 方法 。 文 献 [2] 将 扩展 有 限 元 法 和 扩 
展 等 几何 分 析 分 别 与 粒子 群 优化 和 Jaya 算法 相 结 
合 , 预 测 裂纹 的 位 置 。 文 献 [3] 将 比例 边界 有 限 元 
法 和 机 器 学 习 算 法 相 结合 定量 反 演 结构 中 的 缺陷 。 
作 尖 新 兴 数 值 方法 的 代表 ,数值 流 形 方法 (numerical 
manifold method , NMM) 因 其 特有 的 数学 覆盖 和 物理 
覆 履 系统 已 成 为 求解 裂纹 等 不 连续 问题 的 利器 。 
NNM 求解 裂纹 问题 ( 即 进行 裂纹 问题 正 分 析 ) 的 主 
要 优点 包括 :数学 覆盖 可 不 与 物理 域 的 所 有 边界 
(会 裂纹 ) 一 致 ;可 自然 描述 穿越 裂纹 面 时 位 移 场 的 
不 谋 续 ;通过 恰当 选取 履 盖 函 数 可 准确 表征 裂 尖 区 
物理 场 的 局 部 特性 。 得 益 于 以 上 优点 ,NMM 已 被 广 
泛 岂 用 于 求解 各 类 裂纹 问题 “”。 另 外 ,Piman 神经 
网 络 是 一 种 典型 的 反馈 型 人 工 神经 网 络 , 有 较 强 的 
优化 计算 和 联想 记忆 功能 ,可 以 看 作 是 具有 局 部 记 
忆 导 素 和 局 部 反馈 连接 的 递归 神经 网 络 ""。 CR 
[四 将 Elman 神经 网 络 用 于 缺陷 检测 ,提出 了 一 种 
三 出 结构 的 桥梁 缺陷 自动 检测 与 识别 方法 。 文 献 
[13] 提 出 了 一 种 基于 时 间 序 列 与 Elman 神经 网 络 
的 大 型 复杂 空间 顶 架 结构 损伤 检测 方法 。 

至 今 ,NMM 与 智能 算法 相 结合 已 有 初步 成 果 
(如 文献 [14] 结合 NMM 与 蚁 群 算法 研究 了 边 坡 稳 
定性 ) ,但 其 在 裂纹 反 演 分 析 方面 的 工作 尚未 见报 
道 ,本 研究 将 NMM 与 Elman 神经 网 络 相 结 合用 于 裂 
纹 识别 ,在 探讨 方法 可 行 性 和 精度 的 基础 上 进一步 
分 析 测 点 布置 .数据 噪声 等 因素 对 网 络 预测 能 力 的 
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1 裂纹 反 演 模型 简介 


裂纹 反 演 是 一 类 典型 的 反问 题 。 这 类 问题 是 在 


结构 某 些 关键 点 ( 即 测 点 ) 的 响应 (如 人 位移、 应力 和 
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应 变 等 ) 已 知 的 情况 下 ,通过 反 演 算法 获得 与 该 结构 
有 关 的 一 些 未 知 的 物理 或 几何 特性 。 

为 探测 结构 中 的 裂纹 几何 信息 ,需要 结合 相关 
优化 算法 寻找 合适 的 参数 模型 估计 结构 关键 点 的 响 
应 ,并 与 结构 真实 响应 进行 比较 ,使 得 二 者 之 间 误 差 
最 小 化 , 见 图 1。 

参数 模型 可 通过 下 组 向 量 表示 为 

a-71o,,0,,"",a,)eR" (1) 
Hep n 为 待 反 演 参 数 个 数 。 


Y 


图 1 含 裂纹 结构 及 测 点 示意 图 


Fig.1 Schematic diagram of structure 
with crack and the measuring points 
图 1 中 待 反 演 的 裂纹 的 参数 为 
Qua 7 {A TiAl} (2) 
其 中 (X,Y,) 为 裂 类 i 坐标 (i=1,2,3)。 
反 演 分 析 的 目标 函数 为 


ouye DE (3) 
Hp un (o FI u” 分 别 为 结构 关键 点 响应 量 的 估计 
值 与 真实 值 。 


故 对 裂纹 进行 反 演 分 析 , 即 通过 有 关 优 化 方法 ， 
在 给 定 的 结构 体 范围 内 , 获取 最 优 的 参数 癌 
m, pp 
T A Xa Xa Se (4) 
使 得 式 (3) 最 小 化 。 


2 NMM-Elman 神经 网 络 方法 


2.1 数值 流 形 方 法 


2.1.1 NMM 简介 

NMM 最 初 应 用 于 岩 体 中 节理 和 裂隙 问题 求解 。 
该 方法 利用 现代 数学 “ 流 形 " 思想 ( 即将 许多 个 部 分 
重 和 至 的 区 域 组 合 在 一 起 ,去 覆盖 全 部 材料 体 ) ,采取 
2 套 有 限 履 盖 体 系 ( 即 数学 覆盖 和 物理 覆盖 ) 统一 求 


解 连续 与 非 连续 问题 n9 。 数 学 覆盖 由 许多 数学 分 
片 组 成 ,不 必 与 物理 域 的 边界 (包括 外 部 边界 、 内 部 
裂纹 .孔洞 边界 及 材料 界面 等 ) 保 持 一 致 ,只 须 神 盖 
全 部 物理 域 即 可 ;相互 部 分 重 琶 的 数学 分 片 ,在 形状 
选择 时 亦 不 受 约束 ;物理 分 片 是 由 诸多 数学 分 片 与 
物理 域 的 相交 部 分 构成 ,其 集合 则 构成 物理 覆盖 ; 若 
干 个 物理 分 片 的 公共 部 分 就 形成 流 形 单元 (下 文 简 
称 单 元 ) 。 关 于 上 述 基本 概念 的 更 细致 的 介绍 可 参 
阅 文献 [17] 。 
2.1.2. S A X EAR SC 

基于 上 述 概念 ,NMM 可 将 任意 单元 e 的 位 移 场 
近似 函数 表示 为 1 

uw (X) = Xw(X) u, (X) (5) 

sg X = (X,Y); N 为 构成 单元 e 的 物理 分 片 的 个 
EQO 为 各 物理 分 片 忆 所 对 应 的 单位 分 解 权 函 
数 和 生成 尸 的 数学 分 片上 定义 的 权 函 数 相对 应 ; 


u;(X) -P(X) * a, * &(X) - b, (6) 
:a; RI b, 分 别 为 定义 在 已 上 的 常规 未 知 量 列 阵 
FREIE; P(X) HERE, B 
1 0 X 0 Y 0 (0) 
a | © 
0 1 0 X 0 Y 0 -- 


Pp 于 二 维 问题 ,可 表示 为 [1 
| 0 d, 0 o, 0 Ó, 0 M 


(8) 


[ 81, 8, By, Ba] = | /rsin £ ros 2. rsinfsin £, 


Jrsinĝcos £] (9) 


其 中 (r,6) 为 裂 尖 局 部 坐标 系 中 考察 点 的 极 坐标 。 
2.1.3 离散 方程 的 建立 
基于 修正 变 分 原理 , 可 得 到 NMM 求解 线 弹 性 
裂纹 问题 的 总 体 方程 为 ”1 
KU-F (10) 
其 中 :KU IF 2151] 29 ARR EE AB IE LI MATE E 
列 阵 和 等 效 载荷 列 阵 ,具体 形式 见 文献 [17,19] 。 
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2.2 Elman 神经 网 络 


2.2.1 Elman 神经 网 络 简介 

Elman 神经 网 络 由 Jeffrey Locke Elman 教授 于 
1990 年 为 解决 语音 处 理 问题 而 提出 | ,是 一 种 反馈 
型 局 部 递归 神经 网 络 。 该 网 络 与 传统 BP 神经 网 络 
的 网 络 结构 略 有 差异 ,除了 基本 的 输入 层 、 隐 含 层 和 
输出 层 之 外 ,还 在 隐 含 层 中 存在 用 于 构建 局 部 反馈 
的 承接 层 ,其 拓扑 结构 如 图 2 所 示 。 承 接 层 接收 来 
自 隐 含 层 的 反馈 信号 ,每 个 隐 含 层 节 点 都 有 与 之 相 
应 的 关联 节点 。 作 为 一 种 延 时 算 子 ,承接 层 通过 
自身 存储 对 前 一 时 刻 隐 含 层 的 输出 进行 延迟 和 记忆 
并 反馈 给 隐 含 层 。 通 过 这 种 方式 ,Elman 神经 网 络 
对 历史 状态 数据 更 具 人 敏感 性 ,具有 强大 的 处 理 动态 
信息 、 实 现 动 态 反 馈 的 能 力 。 


输入 层 隐 含 层 输出 层 
2 Elman 神经 网 络 结构 图 
Fig.2 Structure of Elman neural network 
2.2.2 Elman 神经 网 络 工作 原理 
基于 图 2, 用 于 描述 Elman 神经 网 络 的 非 线 性 
状态 空间 的 表达 式 为 


y(E) =glwx(k) +B] (11) 
x(k) 2f[w,x,(k) +w,(u(k-1)) +8] (12) 
x(k) 2ox,(k - 1) ex(k -1) (13) 


式 中 :wi 2 e 3 CBE IR] BOR: ws 为 隐 含 层 
与 输入 层 间 的 权重 ;w 为 输出 层 与 隐 含 层 间 的 权 
E; y CR) 为 输出 层 输出 ;xz(5) BIER SS k OG TNT] 
隐 含 层 输 出 ;x.(h) 为 网 络 第 次 迭代 的 承接 层 输出 ; 
u( 上 一 1) 为 第 -1 个 输入 数据 ;B1\B, 分 别 为 输入 层 、 
隐 含 层 偏 置 量 ; a 表示 承接 层 自 反 馈 系 数 ;f( * ) 为 
隐 含 层 的 非 线 性 激活 函数 , 常 采用 Sigmoid PRAEC, 
见 式 (14) ;g( * ) 为 输出 层 的 线性 激活 函数 。 

Ka) Td Ls 


定义 Elman MAE FIARTUH] 77 BOS RO 


(14) 
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E- VD, 0) -y TD 0) 709]. (15) 


HR: E J E KIR SP Y. CR) 29 PI2R k XS 
TAB] HL y CK) HE k GT SE odi h o 

隐 含 层 神 经 元 个 数 岂 由 经 验 公 式 (16) 确 
i? , 即 

h- /m*n*£ (£z1-10) (16) 

HP im n DIAA i o ns GE 为 1 ~10 [H] 
可 调节 和 常数。 选取 时 一 般 采 用 试 错 法 , 即 在 备 选 
值 范 围 内 找 出 使 均 方 误 差 训 最 小 的 h。 


2.3 NMM-Elman 神经 网 络 方法 的 实现 


本 研究 将 NMM 与 Elman 神经 网 络 算法 相 结合 
运用 于 裂纹 反 演 分 析 。 首 先 基 于 NMM 获取 含 裂纹 
缚 烙 在 外 载 作 用 下 相关 测 点 的 位 移 值 ,然后 将 测 点 
E EL Na MARIA 
数据 ,并 通过 逐步 修正 各 层 权 重 等 参数 .调整 网 络 结 

识别 。 该 方法 的 主要 实现 流程 如 图 3 所 示 。 


基于 NMM-Elman 神 经 网 络 


的 裂纹 反 演 识别 


目标 函数 极 小 化 
搭建 Elman 神 经 网 络 模型 


取 让 NMM 的 收敛 性 
与 计算 精度 


布置 测 点 ， 获 取 各 测 点 
"o 的 位 移 响 应 量 


m 
u 


展 训 经 


将 各 裂纹 构 型 下 获取 的 测 点 
位 移 及 对 应 裂 尖 坐标 导入 拟 
构建 的 Elman 神 经 网 络 中 


将 待 预测 裂纹 构 型 下 由 NMM 所 得 的 
测 点 位 移 作为 输入 数据 ， 基 于 搭建 
的 网 络 预测 其 烈 拓 坐标 ， 将 预测 值 
"SERA, IULII RR 


图 3 NMM-Elman 神经 网 络 方法 开展 裂纹 识别 的 流程 图 


Fig.3 Flow chart of crack detection with 


NMM-Elman neural network method 


3 算 例 分 析 


为 验证 NMM-Elman 神经 网 络 协同 方法 的 可 行 
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性 和 精度 ,本 节选 取 2 个 裂纹 识别 问题 (中 心 水 平 裂 
纹 和 星 形 裂纹 ) 进行 探讨 。 计 算 时 统一 取 杨 氏 模 量 
E 21. 0 GPa, 泊 松 比 v=0.3, 且 均 按 平面 应 变 状态 
考虑 。 
3.1 ERRAR Hbk ERARA 

如 图 4 FER, EJE R a 2h 28.0 m, 宽 2b = 
3.0 m, 板 内 含有 一 条 长 度 为 2a 的 中 心 水 平 裂 纹 。 
板 项 部 受到 o =1.0 MPa 的 拉 应 力作 用 , 板 底 边 的 竖 
向 位 移 和 底 边 左 端 的 水 平 位 移 均 为 零 。 


到 4 单 向 拉 伸 作用 下 含 中 心 水 平 有 裂纹 的 矩形 板 


Fig.4 A rectangular plate with a central horizontal 


crack under uniaxial tension 

3.1.1 NMM 精度 验证 

对 线 弹 性 断裂 力学 问题 ,一 般 通 过 校 核 断裂 参 
数 的 精度 来 验证 方法 的 准确 性 。 据 此 ,此 处 重点 考 
察 应 力 强度 因子 的 NMM 仿真 精度 。 

IRR B 为 例 , 根 据 文献 [23] , 当 h/b zm 1.5 
时 , 裂 尖 应 力 强度 因子 Ki 的 计算 公式 为 

K, = Fovma (17) 

其 中 环 为 有 限 宽 度 修 正 系 数 ,本 算 例 中 取 呈 ) 


1 | +0. 326 (4) 


dis i 。 (18) 
by 
NMM 模拟 时 ,采用 交互 积分 方法 "求解 裂 尖 
应 力 强度 因子 。 
1) 收敛 性 测试 


HL a =0.5m, 分 别 采用 3 种 由 正方 形 数学 单元 
(数学 单元 边 长 ! 分 别 为 0.25 m,0. 15 m,0. 125 m) 
组 成 的 数学 覆盖 对 板 进行 离散 , 剂 分 后 得 到 的 单元 
数量 分 别 为 432、1 086 和 1 632 ,相应 的 离散 域 见 图 
5,NMM 得 到 的 K, 及 对 应 的 参考 解 列 于 表 1。 可 以 
发 现 随 着 覆盖 的 细 化 ,NMM 解 逐渐 收敛 于 参考 解 。 


(a) [-025 m (b) /-0.15 m 


图 5 NMM 离散 域 
Fig.5 NMM discretization 


(c) [-0.125 m 


1 


表 1 不同 覆 盖 尺 寸 下 的 五 (单位 :MPa . mz ) 


Tab.1 K, at different cover size ( unit; MPa * m?) 


NMM 解 
T 参考 解 
2432 单元 1 086 单元 1632 单元 
Q 3234 1.3317 1.3339 1.3344 


02) 裂纹 长 度 取 不 同 值 时 的 计算 精度 
分 别 取 a=0.3,0.4,0.5,…,1.0 m, 并 统一 采用 
l Eg 125 m 对 应 的 数学 覆盖 进行 分 析 。 图 6 给 出 了 
dip NMM 求解 的 K 结果 及 相应 的 参考 解 。 可 以 看 


PNMM 仿真 值 与 参考 解 非 常 吻合 ,表明 NMM 模拟 


1 


- m?) 


J (MPa 


K, 


03 04 05 


0.6 0.7 0.8 0.9 1.0 
裂纹 半 长 wm 


图 6 不 同 裂 纹 长 度 时 的 应 力 强度 因子 

Fig.6 Stress intensity factors at various crack lengths 
3.1.2 Elman 神经 网 络 的 训练 和 预测 

1) 网 络 的 搭建 和 训练 

在 对 Elman 神经 网 络 开展 训练 时 , 将 各 裂纹 长 
度 下 对 应 的 测 点 位 移 值 作为 输入 数据 ,将 裂 尖 4A、B 
的 坐标 ( 即 X, Y, XI Y) 作为 输出 数据 。 其 中 ， 
测 点 的 位 移 通 过 NMM 正 向 分 析 获 得 。 据 此 ,采用 
120.125 m 的 数学 覆盖 先后 对 c =0.3,0.4,0.5,…， 
1.2m 等 10 种 情况 进行 模拟 。 初 步 按 图 7 的 方式 布 
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置 14 个 测 点 ,各 测 点 的 信息 列 于 表 2。 根 据 2.2 节 ， 
Elman 神经 网 络 的 输入 层 和 输出 层 的 神经 元 数 可 确 
定 为 14 个 和 2 个 ,再 根据 式 (16) 并 经 过 试 算 , 取 隐 
含 层 神经 元 数量 为 4, 最 大 训练 次 数 为 1 000 次 , 训 
练 的 目标 均 方 误差 设置 为 0.000 000 1, 损 失 函 数 随 
训练 次 数 的 迭代 收敛 曲线 见 图 8, 可 发 现在 332 次 训 
练 后 满足 误差 要 求 。 为 判断 网 络 的 学 习 效 果 , 首 先 
采用 已 有 训练 数据 对 网 络 进行 测试 ,相关 结果 列 于 
表 3 和 表 4, 可 发 现 预测 值 与 真实 值 非常 接近 ,表明 
搭建 的 网 络 具有 良好 的 自学 习 能 

R2 测 点 坐标 


Tab.2 Coordinates of measuring points 


ASO XAR 了 坐标 | 编号 Xads Y 
1 0.000 00 1.333 33 8 2.000 00 8.000 00 
2 0.000 00 2.666 66 9 3.000 00 8.000 00 
3 0.000 00 4.000 00 10 3.000 00 6.666 66 
4 0.000 00 5.333 33 11 3.000 00 5.333 33 
5 0.000 00 6.666 66 12 3.000 00 4.000 00 
6 0.000 00 8.000 00 13 3.000 00 2.666 66 
7 1.000 00 8.000 00 14 3.000 00 1.333 33 


图 7 测 点 的 布置 


Fig.7 Layout of measuring points 


107 : 
最 优 一 
un 4 
Ie 训练 输出 | 
(332,9.917 4X 107) : 
10* | 
10* 


0 50 100 150 200 250 300 350 
训练 次 数 


到 8 ”训练 迭代 收敛 曲线 


Fig.8 Training iteration convergence curve 
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表 3 Elman 神经 网 络 对 裂 尖 4 坐标 的 训练 输出 结果 


Tab.3 Training results of Elman neural network for crack tip 4 


裂纹 半 长 目标 训练 输出 目标 训练 输出 
a/m X, X, Y, Y, 
0.3 1.2000 — 1.1995 — 4.0000 4.000 0 
0.4 1.1000 — 1.1012 — 4.0000 4.000 0 
0.5 1.0000 | 0.9996 — 4.0000 4.000 1 
0.6 0.9000 | 0.8992 4.0000 3.999 9 
0.7 0.8000 0.7999 4.0000 4.000 1 
0.8 0.7000 0.7010 4.0000 4.000 0 
0.9 0.6000 0.5998 4.0000 3.999 9 
1.0 0.5000 | 0.4997 4.0000 4.000 0 
1.1 0.4000 — 0.4001 4.0000 4.000 0 
pwE 0.3000 | 0.3001 4.0000 4.000 0 


ÈA Elma 神经 网 络 对 裂 尖 B 坐标 的 训练 输出 结果 


TaB 汉 Training results of Elman neural network for crack tip B 


训练 输出 目标 训练 输出 


Ga 长 — Bi 


a/m X, X; Yp Y, 

i 4.3 1.8000 — 1.8005 — 4.0000 4.000 0 
Sen 1.9000 — 1.89088 — 4.0000 4.0000 
c *0.5 2.0000 | 2.0004 4.0000 3.999 9 
C 70.6 2.1000 — 2.1007 — 4.0000 4.000 1 
0.7 2.2000 — 2.2000 — 4.0000 3.999 9 
2S. 8 2.3000 | 2.2990 — 4.0000 4.000 0 
i 30.9 2.4000 2.4002 4.0000 4.000 1 
—1.0 2.5000 | 2.5003 — 4.0000 3.999 9 
Qu. 2.6000 | 2.5999 4.0000 4.000 0 
15 2.7000 2.6999 4.0000 4.000 0 


2) 测 点 布置 对 预测 结果 的 影响 

测 点 布置 方案 :方案 1 ,布置 如 图 7 所 示 14 个 测 
点 ;方案 2 ,布置 除 6.9 外 12 个 测 点 ;方案 3 ,布置 除 
1,14 外 12 个 测 点 ;方案 4, 布 置 除 3.6.9、12 外 10 
个 测 点 ;方案 5 ,布置 除 1.2、13、14 外 10 个 测 点 ; 方 
案 6 ,布置 除 1.3.6.9、12、14 外 8 个 测 点 。 

并 取 a=0.3,0.4,0.5,…,1.2 m 等 10 种 情形 
下 测 点 的 位 移 作 为 输入 数据 对 网 络 进行 训练 ,在 此 
基础 上 ,将 4a=1.3 m 时 的 测 点 位 移 值 作为 网 络 的 输 
入 数据 以 预测 裂 尖 ALB 的 位 置 。 为 合理 确定 各 测 
点 布置 方案 下 隐 含 层 神经 元 的 个 数 ,在 搭建 Elman 
神经 网 络 时 ,进行 了 反复 试 算 , 结 果 表 明 , 当 隐 含 层 
神经 元 个 数 均 取 4 时 ,各 方案 下 的 预测 误差 相对 最 


低 。 相 关 的 预测 结果 见 表 5 ~ 表 6。 可 以 看 出 , 测 点 
的 布置 对 YXs 和 Ys 的 预测 精度 影响 很 小 ,对 XX 的 
影响 则 较 明显 ,这 可 能 与 裂 尖 4 距离 板 左 下 角 点 的 
距离 较 近 而 该 点 施加 了 XX 向 的 位 移 约束 有 关 。 进 一 
步 观察 表 5 和 表 6 还 可 发 现 , 在 6 种 布置 方案 中 , 方 
案 1 的 误差 最 小 。 

表 5 各 测 点 方案 下 裂 尖 4 的 坐标 预测 输出 


Tab.5 Coordinate prediction output of crack tip 


A under different measuring-point layouts 


目标 预测 相对 目标 预测 相对 
X, X, 误差 /% Y, Y, 误差 /% 


不 


1 0.2000 0.1960 2.000 4.0000 4.0000 0.000 


2 0.1953 2.350 3.996 9 0.077 
3 0.1889 5.550 3.999 5 0.013 
4 0.2049 2.450 4.0000 0.000 
5 0.1911 4.450 3.999 1 0.023 
6 0.2046 2.300 4.0002 0.005 


表 6 各 测 点 方案 下 有 裂 尖 8 的 坐标 预测 输出 
Tab.6 Coordinate prediction output of crack tip B 


under different measuring-point layouts 


目标 预测 相对 Hs 预测 相对 


方案 
si Xp Xp 误差 /% Ys Ys 误差 /% 


1 2.8000 2.8040 0.143 4.0000 4.0002 0.005 


2 2.8053 0.189 4.0027 0.067 
3 2.8111 0.396 4.0004 0.010 
4 2.7951 0.175 3.999 9 0.003 
5 2.8091 0.325 4.0009 0.022 
6 2.7954 0.164 3.999 8 0.005 


3 ) 输 入 数据 噪声 对 预测 精度 的 影响 
为 探索 Elman 神经 网 络 的 容 噪 能 力 , 对 由 NMM 
模拟 获得 的 测 点 位 移 值 添 加 噪声 , 即 “ 
Qua = Puuh | +P) (19) 
其 中 ;gi 为 Elman 神经 网 络 的 输入 数据 ( 即 添加 品 
声 后 的 测 点 位 移 值 ) ;pnw 为 NMM 仿真 得 到 的 测 点 
原始 位 移 值 ;p 为 噪声 强度 。 
网 络 训练 仍 考虑 wa =0.3,0.4,0.5,…,1.2m 等 
10 种 裂纹 构 型 ,预测 时 则 将 a =1.3m 的 含 噪声 的 测 
点 位 移 值 ( 按 前 述 方案 1 布置 测 点 , 隐 含 层 神 经 元 取 
4) 作 为 输入 数据 。 图 9 给 出 了 p 分 别 取 0% 、1% 、 
2% 5% 和 10% 时 各 裂 尖 坐标 预测 的 相对 误差 绝对 
值 随 噪声 值 的 变化 情况 。 可 以 看 出 , 裂 尖 坐标 预测 
值 的 相对 误差 随 噪声 值 的 提升 逐渐 增 大 。 此 外 ,还 


可 发 现 加 对 噪声 最 为 敏感 ,这 可 能 与 裂 尖 4 的 位 置 
及 约束 条 件 有 关 ; 为 外 , 当 噶 声 值 为 1% 时 ,XX 的 相 
对 误差 较 无 噪声 时 更 小 ,这 一 现象 可 能 与 网 络 陷 人 
局 部 最 优 有 关 。 当 噪声 值 不 超过 5% 时 ,在 大 部 分 
情况 下 各 坐标 的 预测 相对 误差 均 小 于 1% ,表明 搭 
建 的 Elman 神经 网 络 具 有 很 好 的 容 噪 能 力 。 


一 @ 一 XY 预测 的 相对 误差 绝对 值 
一 会 一 XY, 预 测 的 相对 误差 绝对 值 


e 10 一 TF 一 了 预测 的 相对 误差 绝对 值 
E 一 9 一 7 预测 的 相对 误差 绝对 值 
& 
3 8 
End 
b 6 
E 
Ea 
T- 
> 2 
co 
co H $ : : 
E 噪声 强度 等 级 /% 
e 


图 9 不 同 噪声 下 裂 尖 坐 标 预 测 的 相对 误差 绝对 值 
coy ig.9 Absolute relative errors of the predicted crack 
e tip coordinates under different noises 


co 
302.1 方 板 内 的 星 形 裂 纹 识别 
e 
ON 考察 图 10 所 示 含 星 形 裂纹 的 方 板 。 板 宽 2 下 = 


4. m ,各 分 支 裂纹 长 度 均 为 , 相 邻 分 支 间 的 来 角 皆 
为 加" ,承受 大 小 o - 1. 0 MPa 的 双向 均匀 拉 伸 荷载 
作出, 板 底 边 的 竖 向 位 移 和 底 边 左 端的 水 平 位 移 均 


图 10 双向 拉 伸 荷载 作用 下 含 星 形 裂 纹 的 方 板 
Fig. 10 A square plate with a star-shaped 
crack under bi-axial tension 
3.2.1 NMM 精度 验证 
鉴于 对 称 性 , 仅 考 察 裂 尖 4 和 号 的 应 力 强度 因 
子 值 Ki ,Kis, Kns 的 计算 精度 。 分 别 取 a 20.2W, 
0. 3W,---,0. 6W, {EHK 120. 083 m 的 正方 形 数 
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学 单元 组 成 的 数学 覆盖 对 物理 域 进 行 剖 分 ,NMM 仿 
真得 到 的 应 力 强 度 因子 解 与 对 应 的 参考 解 ”“ A 
表 7, 可 以 看 出 ,NMM 解 与 参考 解 吻合 度 较 高 。 
表 7 RAKE a 取 不 同 值 时 的 应 力 强度 因子 
Tab.7 Stress intensity factors at different crack length a 
K (H/ 裂纹 长 度 a/m 
(MPa - m7) 0.2W 0.3W 0.4W | 0.5W 0.6W 


Kj, (NMM)  Á 0.8567 1.0733 1.2880 1.5183 1.7980 
Ku (参考 解 ) 0.8486 1.0772 1.3087 1.5689 1.8947 
Kjp(NMM) 0.8637 1.1050 1.3632 1.6713 2.0785 
Kip( 人 参考 解 ) 0.8495 1.0824 1.3261 1.6107 1.9769 
Kus(NMM) 0.0006 0.0019 0.0076 0.0237 0.0577 
Kus (参考 解 ) 0.0005 0.0031 0.0111 0.0298 0.0656 


3.2.2 Elman 神经 网 络 的 搭建 和 训练 

1) 网 络 的 训练 

搭建 Elman 神经 网 络 时 依次 取 a = 0. 4043 m, 
0. 42/3 m,0. 4443 m,…,0. 58V3 m 等 10 种 情形 , 仍 
然 将 NMM 模拟 得 到 的 测 点 位 移 值 作为 网 络 的 输入 
RGE, RRA B,CD EF 的 坐标 值 作为 输出 数据 。 
在 NMM 仿真 时 均 取 1=0.083 m( 图 11 给 出 了 a= 
0. 50/3 m 的 离散 域 ,单元 数 为 2476) ,布置 12 个 测 
点 如 图 11 所 示 。 由 此 ,网 络 的 输入 层 、 隐 含 层 和 输 
出 层 的 神经 元 个 数 分 别 为 12 .7( 由 试 算得 出 ) 和 6， 
设置 最 大 训练 次 数 为 1000 次 ,训练 的 目标 均 方 误差 
为 0.000 000 1 ,损失 函数 随 训 练 次 数 的 迭代 收敛 曲 
线 见 图 12。 可 知 ,在 迭代 299 次 后 训练 满足 误差 要 
求 。 随 后 采用 训练 数据 对 网 络 进 行 测试 , 表 8 ~ 表 9 
4p 9 d Y 24e ALB 的 预测 结果 ,表明 搭建 的 
Elman 神 经 网 络 具 有 很 好 的 自学 习 特 性 。 


(0.84.0) (1.64.0) (2.4 4.0) (3:2 4.0) 


图 11 NMM 离散 域 及 测 点 布置 (a =0. 50/3 m) 
Fig.11 NMM discretization and measuring 


points layout when a = 0. 50/3 m 
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训练 迭代 收敛 曲线 
Fig.12 Training iteration convergence curve 
表 8 Elman 神经 网 络 对 裂 尖 4 坐标 的 训练 输出 结果 


Tab.8 Training results of Elman neural 


network for crack tip A 


作为 输入 值 以 预测 裂 尖 坐标 ,相应 的 结果 列 于 表 10 
和 表 11。 可 以 看 出 ,对 所 有 裂 尖 , 预测 的 相对 误差 
最 大 值 均 不 超过 0.2% ,表明 搭建 的 Elman 神经 网 


络 具 有 较 高 的 泛 化 能 力 。 


表 10 Elman 神经 网 络 对 有 裂 尖 坐 标的 
预测 值 和 相对 误差 (a =0. 60/3 m) 


Tab. 10 Prediction and relative error of Elman neural 


network for crack-tip coordinates (a =0. 60/3 m) 


裂纹 长 度 om Hs X, 训练 输出 XY。 HERY, WRH Ya 
E 0. 40/3 2.6928 2.692 8 2.0000 1.999 9 
| "n. 42/5. 2.7215 2.7214 2.0000 2.0001 
Cf aA 27621 — 2.7238 — 2.0000 — 2.0000 
e 9. 46/3 2.796 7 2.796 8 2.0000 2.0000 
t Y». 48/3. 2.8314 2.8312 2.0000 1.999 9 
CO so 2.8660 2.8663 2.0000 1.9999 
^ à. 52/5. 2.9007 2.9004 2.0000 2.0001 
- 54/3 2.9353 2.935 4 2.0000 1.999 9 
x H 56/3 2.969 9 2.9700 2.0000 2.0000 
a, 58/3 3.004 6 3.004 6 2.0000 2.0000 
( 表 9 Elman 神经 网 络 对 裂 尖 B 坐标 的 训练 输出 结果 
C : Tab.9 Training results of Elman neural 
( network for crack tip B 
AAKE a/m 目标 Xs 训练 输出 Xs 目标 Y。 训练 输出 Ys 
0. 40/3 2.3464 2.3464 2.6000 2.6000 
0. 42/3. 2.3637 2.3637 2.6300 2.6299 
0. 44/3. 2.3811 2.3811 2.6600 2.6601 
0. 46/3. 2.3984 2.3983 2.6900 2.6900 
0. 48/3. 2.4157 2.4156 2.7200 2.7198 
0. 50/3. 2.4330 2.4331 2.7500 2.7503 
0. 52/3. 2.4503 2.4503 2.7800 2.7798 
0. 54⁄3 2.4677 2.4677 2.8100 2.8101 
0. 56/3 2.4850 2.4849 2.8400 2.8399 
0. 58/3 2.5023 2.5023 2.8700 2.8700 
2) 网 络 的 预测 精度 


分 别 取 a =0. 60/3 m 和 a 20.62/3 m ,基于 已 拱 
建 完 成 的 Elman 神经 网 络 , 将 NMM 得 到 的 测 点 位 移 


RR 坐标 目标 值 预测 值 相对 误差 /% 
A X, 3.0392 3.0415 0.076 
Y, 2.0000 2.0002 0.011 
B Xp 2.5196 2.5204 0.032 
Y, 2.9000 2.9007 0.023 
C Xe 1.480 4 1.4798 0.043 
Ye 2.900 0 2.900 2 0. 008 
D Xp 0.9608 0.9600 0. 080 
Yp 2.0000 1.9998 0.010 
E Xe 1.480 4 1.479 4 0.065 
Y; 1.1000 1.0994 0.059 
F Xr 2.5196 2.5205 0.037 
Yp 1.1000 1.099 5 0.044 
表 11 Elman 18S [9] £& 9E 23 Ze 5E B] 
预测 值 和 相对 误差 (a =0. 62/3 m) 
Tab.11 Prediction value and relative error of Elman neural 
network for crack-tip coordinates (a 20. 62/3 m) 
RR 坐标 目标 值 预测 值 相对 误差 /% 
4 X4 3.0739 3.0756 0.056 
Y, 2.0000 1.999 7 0.015 
B Xp 2.5369 2.5372 0.012 
Ys 2.9300 2.928 7 0.046 
C Xe 1.463 1 1.4625 0.036 
Ye 2.9300 2.9281 0.066 
D Xp 0.926 1 0.9264 0.031 
Yp 2.0000 1.999 6 0.018 
E Xe 1.463 1 1.4626 0.031 
Y, 1.0700 1.0715 0.138 
F Xp 2.5369 2.5376 0.026 
Yp 1.0700 1.0712 0.115 


进一步 探讨 输入 数据 噪声 对 预测 结果 的 影响 。 


图 13 ~ 图 14 分 别 给 出 了 噪声 值 依 次 为 0% 、1% 、 


2% 5% 和 1096 时 ,Elman 神经 网 络 对 真实 裂纹 a 取 
0. 6043 m 和 0. 62/3 m 的 反 演 结果 。 不 难 发 现 ,网 络 
预测 的 误差 随 噪声 强度 等 级 的 提升 而 增 大 ;此 外 ,对 
所 取 的 噪声 强度 ,预测 相对 误差 绝对 值 保 持 在 8% 
以 内 , 且 绝 大 部 分 低 于 2% ,说 明 网 络 容 噪 性 能 好 ， 
稳定 性 强 , 个 别 裂 尖 误差 稍 大 (如 X, 和 Yi) 的 结果 
可 能 与 裂纹 的 位 置 及 约束 情况 有 关 。 整 体 而 言 , 搭 


建 的 网 络 模型 可 以 满足 预测 要 求 。 
—0- a-0.604/3 m， 马 预测 的 相对 误差 绝对 值 
A q=0.62V3 m， 了 XY 预测 的 相对 误差 绝对 值 
—V— a-0.604/3 m， 了 ,预测 的 相对 误差 绝对 值 
8 r-e- a-0.6243 m，Y, 预 测 的 相对 误差 绝对 值 
一 4 a-0.604/3 m，XX. 预 测 的 相对 误差 绝对 值 
7 | P a=0.62V3 m，X 预 测 的 相对 误差 绝对 值 
一 全 a-0.604/3 mm， 总 预测 的 相对 误差 绝对 值 
—k- qa=0.62V3 m， 有 %, 预 测 的 相对 误差 绝对 值 
6 上 会- ao=0.60V3m， 总 预测 的 相对 误差 绝对 值 
—9— o=0.62V3m， 闷 预测 的 相对 误差 绝对 值 
5 L 470.6043 m， 闷 预测 的 相对 误差 绝对 值 
—- ao=0.62V3m， 六 预测 的 相对 误差 绝对 值 
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ig.13 Absolute value of relative error of X-coordinate 
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© prediction under different noises 
c —@- a=0.60V3 m， 了 ,预测 的 相对 误差 绝对 值 
c ê- 4-0.624/3 m. 了 ,预测 的 相对 误差 绝对 值 
8 r5 a-0.6043 m， 了 ,预测 的 相对 误差 绝对 值 
C) |$ o-0.62V3m， 力 预测 的 相对 误差 绝对 值 
7 [4 470.6043 m, Y GU EDGE YR 2 £O] t 
-> a-0.6243 m， 了 .预测 的 相对 误差 绝对 值 
-全 ~ a-0.6043 m， 玉 预测 的 相对 误差 绝对 值 
6 F- a-0.6243 m， 了 ,预测 的 相对 误差 绝对 值 
—$- a-0.6043 m, Y, 对 误差 绝对 值 
5 [-9- 470.6243 m, Y 对 误差 绝对 值 

—— qa=0.60V3 m, Y, 对 误差 

-ou=0.62V3m， 帮 预测 的 相对 误 


相对 误差 的 绝对 值 /% 
ES 
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图 14 不同 噪 声 下 了 坐标 预测 值 的 相对 误差 绝对 值 


Fig. 14 Absolute value of relative error of Y-coordinate 
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4 结 论 


双重 覆盖 系统 的 引入 使 得 数值 流 形 方法 可 高 精 
度 便捷 地 模拟 各 类 裂纹 问题 ,承接 层 的 存在 则 确保 
T Elman 神经 网 络 的 自学 习性 和 泛 化 能 力 。NMM 
可 在 固定 数学 覆盖 上 求解 连续 变化 的 裂纹 问题 的 优 
点 ,很 好 地 迎合 了 Elman 神经 网 络 对 训练 数据 获取 
的 精度 及 效率 的 需求 。 据 此 ,本 研究 发 展 了 用 于 线 
弹性 体内 部 直线 裂纹 反 演 分 析 的 NMM-Elman 神经 
网 络 协同 方法 ,给 出 了 两 者 的 协同 方式 ,通过 2 个 典 
型 算 例 验 证 了 方法 的 精度 和 和 鲁 棒 性 ,并 探讨 了 测 点 
布置 .输入 数据 噪声 等 对 预测 精度 的 影响 ,结果 
明 ,本 研究 的 方法 具有 很 好 的 裂纹 反 演 能 

为 了 便于 与 已 有 文献 中 的 参考 解 对 照 ,本 研究 
的 2 个 算 例 均 选用 了 对 称 裂 纹 , 实 际 上 ,本 研究 的 方 
法 同样 适用 于 非 对 称 直 线 裂 纹 的 反 演 。 同 时 ,本 研 
究 的 分 析 对 象 为 静止 直 裂 纹 ,后 期 将 进一步 围绕 曲 
裂纹 动态 裂纹 等 更 复杂 的 问题 开展 进一步 研究 。 
此 外 ,在 分 析 数 据 噪声 对 Elman 神经 网 络 的 影响 时 ， 
本 研究 对 所 有 测 点 数据 采用 了 一 致 的 噪声 值 , 后 期 
还 将 围绕 更 符合 实际 的 随机 噪声 的 影响 开展 相关 
工作 。 
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